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Abstract 

We consider the d = 1 nonlinear Fokker-Planck-like equation with fractional 
derivatives ^P(x,t) = D-§^ [P (x,t)] p '. Exact time-dependent solutions are 
found for v = j^- (—00 < 7 < 2). By considering the long-distance asymp- 
totic behavior of these solutions, a connection is established, namely q = 
(0 < 7 < 2), with the solutions optimizing the nonextensive entropy char- 
acterized by index q . Interestingly enough, this relation coincides with the 
one already known for Levy-like superdiffusion (i.e., v = 1 and < 7 < 2). 
Finally, for (7,^) = (2,0) we obtain q = 5/3 which differs from the value 
q = 2 corresponding to the 7 = 2 solutions available in the literature {y < 1 
porous medium equation), thus exhibiting nonuniform convergence. 
PACS: 05.60,+w, 05.20.-y, 05.40.+j, 66.10.Cb 
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A great variety of diffusive problems in nature, namely those referred to as normal 
diffusion, are satisfactorily described by the Fokker-Planck linear equation 





dt 



P(x,t) = DV 2 P (x,t) 



(1) 



where P (x, t) is the density of probability in the x = {x%, x 2 , xj} space and D > is the 
diffusion coefficient. Such processes are currently characterized by the fact that (x 2 ) oc t, as 
shown by Einstein in his celebrated 1905 paper on Brownian motion. 

More recently, several works [HI have focused on the same type of linear equation but 
with fractional derivatives. More precisely 



^P (x, t) = D V 7 P (x, t) (-00 < 7 < 2) 



(2) 



where V 7 = J2f=t ~§^r ■ Also, the nonlinear equation with ordinary derivatives has been 0| 



focused on as well. More precisely 



d 



w>-r 



(3) 



(no solutions are known for v < — 1 which are integrable 0]). 

These two generalized Fokker-Planck equations have been used to study anomalous Levy- 
like diffusion as well as correlated-like diffusive processes in porous media [PHT7|| . The present 
paper addresses the unification of both equations as follows: 



^P(x,t) = DV~<[P(x,t)Y (-oo< 7 <2) 



(4) 



We will restrict ourselves to the d = 1 case. More specifically, we are interested in normalized 
scaled solutions of the type 

1 . 



P(x,t) 



0(t) 



X 



0(t) 



(5) 



Inserting this form into Eq.(01) (and, without loss of generality, setting D = 1) we obtain: 



0(t) 2 



±F{z) + zF{z) 



(P 



<f> (t) v+1 dz^ 



[F(z) 



(6) 



where we have used the generic property 

^F(ax) = a s ^F(z) (Sell) (7) 

with z = ax. This basic property holds not only for the ordinary derivative but also for all 
fractional operators we are aware of. By choosing the ansatz 



where k is an arbitrary constant, we obtain 



k (8) 



(t) = (9) 

(kit + k 2 ) 

£- lF{z)r = k jL lzF{z)] (io) 

with k\ = — (7 + v — 1) k, k 2 being another arbitrary constant. Finally making an integra- 
tion we obtain 



-{F{z)] u = kzF{z)+c (11) 



where c is another arbitrary constant. 

Thus far it has not been necessary to specify the fractional operator we refer to. Indeed, 
several fractional generalizations exist for the ordinary derivative, namely the Riemann- 
Liouville |IJ (based on Laplace transform), Weyl ffl] (based on Fourier transform) and Caputo 
|T8| (also based on Laplace transform) ones. From now on we will use the Riemann-Liouville 
operator, since it is for this one that it has been possible to find new exact solutions. In this 
case we will work with the positive x axis and, later on, we will use symmetry to extend the 
results to the entire real axis (we are working, in other words, with 507). Also, we will use 
the following generic result |T3 (see Appendix): 



Dt \x a (a + bxf] = Q * r t a + 1 ] x *-s < a + hx fs (12 ) 
x l K i r [a + 1 - 5} K ' v ' 

with D 5 X = and 5 = a + (3 + 1. By defining g(x) = x» (a + bx)~ and A = a (l — ~ ) — 5, 
and rearranging the indices, Eq . (p!2|) can be rewritten as follows: 



Dl [g{x)\ l 



V [a + 1] 



a s x x g (x) . 



r [a + 1 - 5] 

Using this property in Eq. (|TT|) and, for simplicity, choosing c = 0, we find 

(2-7)7 



a 



1 - 2 7 
7 2 - 3 7 + 2 
1 - 2 7 



2-7 
^ = . 

1+7 

These results allow us to write the solution in the form 

a r z^+v 



P(x,t) 



7+1 



.4 



(Ihlt)^^ V{l + bzf 7 

r(/3) 



T(a + 1)) 



1 + 7 
1-27 



X 



(13) 

(14) 
(15) 
(16) 

(17) 

(18) 
(19) 



where h is an arbitrary constant (to be taken, later on, as ±1 according to the specific 
solutions that are studied) and where, without loss of generality, we have set k<i = and 
a = 1. Indeed, the ki constant can be incorporated into a shift of the origin of time, and a 
can be incorporated into the normalization constant A. We also mention the exact solution 

F{z)ocz^ (20) 

that is not normalizable. Several regions will have to be considered, namely 

(21) 



1 1 

-oo<7<-1,-1<7<0,0<7<-,-<7<1,1<7<2 



We start with the region — oo < 7 < — 1 for which, again without loss of generality, we 
can choose b = — 1. The normalization condition implies 



1 r 



.4 



^7(7+1) 



'X-z 



1_ 7 2 



1 

I-27 



dz = 2A ■ 



r 


V-7+1" 


r 


r 7(7-2)] 




1-27 


27-1 



1 



(22) 



r[i- 7 ] 

See Fig. 1. Also, as we can see from the limits of the different regions (pl|), we will have to 
consider different particular cases, namely 7 = —1,7 = 0,7 = 1/2,7 = 1 an d 7 = 2. 



FIGURES 




7+1 7+1 

FIG. 1. (\ki\t)~t' 2 -~t +1 P(x,t) versus z = x/{\k\\t)~< 2 -^+ 1 for 7 < 0; 7 = _ corresponds 
to a distribution everywhere vanishing, except at the abcissa being ±1 where it diverges. For 
7 < —1 (—1 < 7 < 0) the distribution vanishes (diverges) at z = dbl. 

Let us start with 7 = — 1 and arbitrary v. The corresponding equation is 



d_ 

Of 



X 



P(x,t)= J [P(y,t)Ydy (23) 




To solve it let us go back to Eq. (|T0l); after derivation with respect to z, we obtain 



k^^- = [F(z)r (24) 

We are not going to treat this equation in detail; we rather limit ourselves to remark that 
the value 7 = —1 corresponds to v — > ±00 in the curve (16). In the region — 1 < 7 < the 
probability density has a compact support, on the edges of which it diverges. See Fig. 1. 

Let us now address the < 7 < 1/2 region (where 6=1). In the limit 7 = (7 = 
1/2) integrability fails at infinity (at the origin); no such problems exist for < 7 < 1/2. 
Normalization implies 



It is easy to show that P (x, t) achieves a maximum (see Fig. 2) at 

7 

Z (1 - 2 7 ) 

Let us now address the 7 = limit. It corresponds to the equation 



(26) 



^P(x,t) = [P(x,t)Y (27) 

that can be resolved analytically for arbitrary v. To obtain this solution it is convenient to 
go back to Eq. (p~Q|) . It follows that 

B 



F(z) 



1 + cz x ~ v 



(28) 



Therefore, in the z — > oo (or, equivalently x — > oo) limit, we have that F(z) ocl/zifi/>l 
and -F(z) is a constant if v < 1. The u = 1 case needs specific discussion and we obtain 
F(z) oc zfe -1 . It is worthy reminding that the 7 = solutions cannot be considered as 
distributions of probabilities since they are not normalizable. 

The 7 = 1/2 limiting case corresponds to the following linear equation: 

lp( x ,t) = -j^P(x,t) (29) 

This equation can be solved by using the Laplace transform on both t and x. It can also 
be solved by taking the limit 7 — > 1/2. We have followed this procedure and, after tedious 
though straightforward calculations, we obtain k ~ |(1 — 27) and A ~ [^(1 — 27)] _1 / 2 /2, 
which, replaced into Eq. (|17D, yield 



P 1/2 (x, t )^limP 7 (x, t )= 4 ^ t2 ^ 2) ^ , (30) 

which is a distribution of the Poisson type. Let us stress that the above distribution can 
indistinctively be obtained by solving the fractional differential equation for 7 = 1/2, or by 
taking the 7 > 1/2 and the 7 < 1/2 solutions and then considering the 7 — > 1/2 limit. 
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Let us now focus on the 1/2 < 7 < 2 region (where b = — 1). The solutions strictly 
vanish inside an interval which contains the origin. Outside this interval, the solutions are 
everywhere finite if 1/2 < 7 < 1, whereas they diverge if 1 < 7 < 2 (see Fig. 2). 




7+1 7+1 
FIG. 2. (|fci|i) t 2 -7+i p (x, t) versus z = x/(\k\\t) t 2 — y+i for < 7 < 2; 7 = + corresponds to a 

distribution which vanishes everywhere. For < 7 < 1/2 the distributions are defined in the entire 

real abcissa axis; for 1/2 < 7 < 1 the distributions vanish within the (—1, 1) abcissa interval; for 

7 > 1, a divergence exists at the abcissa ±1 (vertical dashed asymptotes). 



The solutions corresponding to the 1/2 < 7 < 1 region are as follows: 



P(x,t) 



A 



7+1 

(H)t 2 -7+1 



-l + z 



1-7 2 



1-27 



(31) 



Normalization implies 
2A 



00 r .7(7+1) 



1 



[-l + z 



,l-7 2 



1-27 



dz = 2A ■ 



r[ 7 ]r 


7(7 
. 2 7 


-2)1 
-1 


r 


[7(7+1)" 
27-1 





from which A is uniquely determined; finally, k is obtained from A by using Eq. 
Fig. 2. 



(32) 
See 
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Let us now focus the special case 7 = 1. In this case the equation becomes 

l^'^i" <*■«>"■ < 33 > 

Its generic solution of the form indicated in Eq. (|5])) is: 

[F(z)] u = kzF(x,t) + c, (34) 
which implicitly determines F(z). The solution corresponding to c = is 

F (z) oc . (35) 

In the region 1 < 7 < 2 we have the same analytic solution that we had in the region 
1/2 < 7 < 1 (i.e., Eq. ([H])); however, at the point z — 1, a divergence is now present 
(see Fig. 2). It is clear that this solution cannot be used without appropriate asymptotic 
considerations for 7 = 2, since for 7 — > 2, v — > 0. 

Let us now specifically focus on the possible x — ► 00 tail of P(x,t) for arbitrary t. For 
— 00 < 7 < the support is compact, hence there is no tail. For < 7 < 2, we obtain, using 
either Eq. flT7|) (for < 7 < 1/2) or Eq. (J3T]) (for 1/2 < 7 < 2), the following asymptotic 
behavior: 

7(7+1) 

1 7(7 + l) + 7 2 -l t7 2 "7+l 

P (x, t) ~ 7+1 ^ ^ ~ (36) 

^7 2 -7+l X 



We can easily verify that the exponent [7 (7 + l)]/[7 2 — 7 + 1] monotonically increases from 
zero to + 1 when 7 increases from zero to 1+ 2 V T Also, we verify that, for < 7 < 1, both 
{\x\) and (x 2 ) diverge; for 1 < 7 < 2, only the latter does. For 7 < 0, all momenta are finite 
since the support is compact. Finally, in all cases, (x) vanishes because of symmetry. 

Let us now address the solutions for v = 1 and arbitrary 7. The equation to be solved 
becomes 

§i p <■*•*> = £ p <■*•<> < 37 > 



This equation can be solved using Laplace transform and the solutions are discussed in P0 



Moreover, we can see (in Fig. 3), that the point (7, v) = (1/2, 1) is at the crossing of two 



solvable lines. Since solution (30) has been found as the limit of either one of these lines, 
it seems reasonable to conjecture that the same solution is found as a limit along any curve 
through that point. 



v 



1 3r 














-1 

]_ 


1 2 



Y 



FIG. 3. Curves in the (7,^) space on which exact solutions are now available: v = 1 [1,20], 
7 = 1 (present work), 7 = 2 [3], v = (present work). The horizontal dashed line corresponds 
to the 7 — ► —00 asymptote; the vertical dashed line corresponds to the v — > ±00 asymptote. 



Let us summarize the present work. We have addressed a generic Fokker- Planck- like 
diffusive equation, namely the one-dimensional case of Eq. (^), and have looked for exact 
scaled solutions of the type in Eq. (|3p. In the (7, v) parameter space, the solutions corre- 
sponding to 7 = 2 and v > — 1 (porous medium equation) as well as to < 7 < 2 and v — 1 
are available in the literature, as already mentioned. We are now exhibiting exact solutions 
along two new lines, namely the line 7 = 1 and arbitrary v, and the line indicated in Eq. 
(16) (i.e., v — (2 — 7)/(l + 7) with —00 < 7 < 2). For the latter, we observe on Eq. (36) 
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that the spatial asymptotic behavior is characterized by the exponent 1 + 7, which exactly 
coincides with that corresponding to Levy superdiffusion. This is a remarkable result, since 
the present solutions concern a nonlinear Riemann-Liouville- fractional differential equation, 
! and not the usual linear Fourier-fractional one, whose solutions are known to essentially 
be Levy distributions. It would no doubt be interesting to know whether the same behavior 
is obtained no matter the value of v (0 < 7 < 2). 

Let us finally mention a connection between the present problem and the solutions 
obtained from the optimization, under appropriate constraints (normalization and finite 
g-expectation value of x 2 in the interval (—00,00)), of the nonextensive entropy |2l|j22]| 
S q = [1 — / dx p{x) q ]/[q — 1]. It has been shown that these optimizing distributions precisely 
coincide with the solutions of the present diffusive problem for 7 = 2. It comes out that 
q = 2 — v (y > —1) 0. Along the line indicated in Eq. (16), the exact solutions of the en- 
tropic optimization problem and the present diffusive one do not coincide for arbitrary value 
of x. However, comparison of the x — > 00 asymptotic behaviors is possible [p^fl . Indeed, by 
identifying the behavior exhibited in ( |3T)| ) with the behavior l/\x\ 2 ^ q ~~ 1 ' > obtained for 
the entropic problem, we obtain 

? = I XT (°<^ <2 )' ( 38 ) 
7 + 1 

which, as commented above, precisely reproduces the connection established for Levy distri- 
butions [[/]]. By using Eq. (16), this relation can be rewritten as follows: 

5 + 2v , . , , 

q = — ^— (0<u<2). (39) 

(We remind that the distributions for v > 2 and v < —1, i.e., 7 < 0, have compact support). 
The present nontrivial solution provides, for (7,//) = (2,0), q = 5/3, whereas the porous 
medium equation solution q = 2 — v provides q = 2. This discrepancy exhibits that the 
point (7, v) = (2, 0) is a singular one, at least within the fractional derivative that we have 
adopted in this work. 

Another point worth to be mentioned is that we have compared the present solutions 
with those optimizing S q defined in the interval (—00, 00) and using finite g-expectation 

10 



for x 2 (whereas the g-expectation value of x vanishes). This is appropriate since, through 
symmetry, we have extended the solutions that we have found in the interval (0, oo) to the 
entire real axis. Another possibility would of course be to compare the present results in the 
positive real semi-axis with those optimizing the entropy S q defined in the same semi-axis 
and using finite g-expectation value of x. If we did so, relation (39) would be replaced by 
?=(7 + 2)/(7 + l) (see also @||). 

Finally, last but not least, now that we have finished the presentation of the various 
exact solutions that emerged within the present work, it is worthy to mention that it would 
be very welcome the discussion of the stability of such solutions. More precisely, if we start 
at t = with an arbitrary distribution P(x, 0) and make it evolve through the present 
differential equation, what would be the t — > oo asymptotic distribution P a (x, t)l (a stands 
for attractor, in the space of the distributions). For instance, if the evolution is determined 
through the convolution product (i.e., a linear fractional-derivative Fokker-Planck equation, 
using a Fourier-based definition of fractional derivative), then the standard and the Levy- 
Gnedenko central limit theorems apply, and consequently the attractor P a (x,t) is either a 
Gaussian or a Levy distribution (respectively when the second cumulant is finite or infinite). 
If the evolution is instead det! ermined through a nonlinear integer-derivative Fokker-Planck 
equation like the one considered in fU, then P a (x, t) is given (as numerical verifications have 
shown) by the distributions which optimize the nonextensive entropy, where x scales with a 
simple function of t. If the time evolution is obtained, as sometimes done, through recursive 
use of maps p5| , P a (x,t) can present a variety of shapes depending on the specific map 
which is used. Finally, in our present case (nonlinear fractional-derivative Fokker-Planck 
equation using a Laplace-based definition of fractional derivative), the solutions we have 
found might well be P a (x,t). This point, however, deserves analysis on its own. 



APPENDIX 

We give a short review of the property of the fractional operator used to solve the 
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nonlinear equation (4). The Riemann - Liouville operator is used in many applications of 
fractional calculus. The usual integral representation for this operator is: 

l d n } dtfit) 



D a J (x) 



r (n — a) dx n J (x — t) 



a+l 



(n — 1 < a < n) 



(A-l) 



For the calculations in this paper we have instead used the following equivalent form: 



(A-2) 



r(p + l-a)' 

We have also used the generalized Leibnitz formula for this kind of fractional derivative, 
namely 



D a x [f (x) g (x)} = £(* ) D T n If (*)] D l 19 (x) 



(A-3) 



Let us also mention that, for this operator, the following property holds under Laplace 
transform: 



L[D a J{x)) = s a F{s)-f a -V ( o) 



(A-4) 



for a < 1 and where (0) means fractional derivative calculated in x = 0. 

Let us now show how formula ( |A-3| ) leads to Eq. (|12D, that is used in the text. By 

assuming no restrictions on the parameters a, f3 and 7 and applying the generalized Leibnitz 

rule to the function x a (a + bx) 13 we obtain 

00 / 5 \ 



x [a 



+ bxf =E )D 5 - n [x a )D: (a + bx) 



n=0 



v n 



(A-5) 



After some algebra we obtain 



D° 



x a (a + bxf] = £ 



'5 s 



aT^ jW+l) ( _ ir h n {a + bx) ,~n 



n=0 



,n) T (a + 1 — S + n) 



F(J3 + l-n) 



(A-6) 



A closed form for this series can be achieved if a + (3 + 1 = S. Indeed, by using the Gamma 
function property F (z)T (1 — z) = . n we obtain 



Dt 



x a (a + bxf 



a> , r(a + 1 > 
r (a + 1 - 8) 

12 



a + bx) 



0-5 



(A-7) 



which essentially is Eq. (12). 

For completeness, it is worthy to briefly mention here a recent variation of Riemann- 
Liouville operator that we have mentioned in the text, namely, the Caputo derivative. Its 
definition is: 



(C stands for Caputo). The main advantage with respect to the Riemann-Liouville operator 
is that Caputo derivative of a constant is zero, which is not the case of the Riemann-Liouville 
one. Substantially, this kind of fractional derivative is a formal generalization of the integer 
derivative under Laplace transform. As disadvantage, it exhibits the fact that, whenever 
the derivation index is an integer number, it recovers the usual derivative excepting for an 
additive constant, whereas the Riemann-Liouville operator has no such disagreable property. 

Finally, let us also mention the definition of Weyl fractional derivative. It is based on 
the properties of Fourier transform, and it is defined as follows: 




(A-8) 




E 




(A-9) 



k=- 



oo 



its continuum version being 




1 



/. 



■+oo 



doo (—iu) a f (oo) exp (— iujx) . 



(A-10) 



— oo 



(W stands for Weyl). 



13 



REFERENCES 

[1] K.S. Miller and B. Ross, An introduction to the Fractional Differential Equations (Wiley, 
New York, 1993). 

[2] A.R. Plastino, A. Plastino, Physica A 222, 347, (1995). 

[3] C. Tsallis and D.J. Bukman, Phys. Rev. E 54, R2197, (1996). 

[4] G. Drazer, H.S. Wio and C. Tsallis, Phys. Rev. E 61, 1417 (2000). 

[5] R.Hilfer and L. Anton, Phys.Rev.E, 51, R848, (1995); L. Anton, R. Hilfer, Phys. Rev. 
E, 59, 6819, (1999); R. Hilfer, Fractals 3, 211 (1995); R. Hilfer, in Anomalous Diffusion: 
From Basis to Applications, eds. R. Kutner, A. Pekalski and K. Sznajd-Weron (Springer, 
Berlin, 1999) page 77. 

[6] PA. Alemany and D.H. Zanette, Phys. Rev. E 49, R956 (1994); D.H. Zanette and 
PA. Alemany, Phys. Rev. Lett. 75, 366 (1995); M.O. Caceres and C.E. Budde, Phys. 
Rev. Lett. 77, 2589 (1996); D.H. Zanette and PA. Alemany, Phys. Rev. Lett. 77, 2590 
(1996). 

[7] C. Tsallis, S.V.F Levy, A.M.C. de Souza and R. Maynard, Phys. Rev. Lett. 75, 3589 
(1995) [Erratum: 77, 5442 (1996)]; C. Tsallis, A.M.C. de Souza and R. Maynard, in 
Levy flights and related topics in Physics, eds. M.F. Shlesinger, G.M. Zaslavsky and U. 
Frisch (Springer, Berlin, 1995), p. 269; D. Prato and C. Tsallis, Phys. Rev. E 60, 2398 
(1999). 

[8] A.S. Chaves, Phys. Lett. A 239, 13 (1998). 

[9] M. Buiatti, P. Grigolini and A. Montagnini, Phys. Rev. Lett. 82, 3383 (1999). 
[10] A.I. Saichev, Chaos 7, 753 (1997). 

[11] V. Sechadri and B.J. West, Proc. Natl. Acad. Sci. USA 79, 4501 (1982). 
[12] H. Spohn J. Phys. I 3, 69, (1993). 

14 



[13] O.V. Bychuk and B. O'shugness, Phys. Rev. Lett. 71, 3975, (1993). 

[14] D.E. Strier, D.H. Zanette, and H.S. Wio, Physica A 226, 310, (1996). 

[15] B.J. West, P. Grigolini, R. Metzler and T.F. Nonnenmacher, Phys. Rev. E 55, 99 
(1997); B.J. West and P. Grigolini, in Applications of Fractional Calculus in Physics, 
ed. R. Hilfer (World Scientific, Singapore, 1998). 

[16] P. Grigolini, A. Rocco and B.J. West, Phys. Rev. E 59, 2603 (1999). 

[17] T. Huillet, J. Phys. A 32, 7225 (1999). 

[18] M. Caputo, Elasticitd e dissipazione (Zanichelli, Bologna, 1969); M Caputo and F. 
Mainardi, Riv. Nuovo Cimento (Ser. II) 1, 161 (1971). 

[19] M.Bologna, Derivative of Real Index [in Italian], ed. E.T.S. Pisa (1990). 

[20] R. Metzler and T.F. Nonnenmacher, Phys. Rev. E 57, 6409, (1998). 

[21] C. Tsallis, J. Stat. Phys. 52, 479, (1988). 

[22] C. Tsallis, Braz. J. Phys. 29, 1, (1999) [accessible at 

|http: / /sbt.if. usp.br/ WW W_pages/Journals/BJP/ Vol29/Muml/index.htm[ . 

[23] The distributions obtained through optimization of the nonextensive entropy S q , the 
present ones and Levy distributions do not generically coincide. However, they all ex- 
hibit a power-law at long distances, hence identification of their respective asymptotic 
behaviors becomes possible. 

[24] A.K. Rajagopal and S. Abe, preprint (2000) ||cond-mat/0003304 |. 

[25] M. Ignaccolo and P. Grigolini, A non-extensive approach to the time evolution of Lya- 
punov coefficients, preprint (2000). 



15 



